
* Load data
use "${data_out_folder}/belief_validation.dta", clear

* Generate figures
twoway (hist belief, width(25) lcolor(gs6) fcolor(gs8) percent) (scatteri 0 225 20 225, recast(line) lpattern(solid) lwidth(0.5) lc(maroon)) (scatteri 0 127.6 20 127.6, recast(line) lpattern(dash) lwidth(0.5) lc(gs4)),       xlab(0(50)450, nogrid labsize(medium))  ylab(, nogrid labsize(medium)) xsize(4) ysize(3) xtitle("Belief about the average donation amount ($)", size(medlarge)) ytitle("Percent", size(medlarge)) legend(order(2 "Average donation" 3 "Average belief about the donation") pos(12) row(1) size(medium))

graph export "${out_folder}/figures/fig_A5.pdf", replace

* For the paper:
gen underestimator = belief <= 225
summ underestimator
gen wedge = belief - 225
summ wedge
